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ABSTRACT 

We present an example of using a combined ray trace and diffraction modeling code to simulate effects of 
objective-lens tilt in an optical data storage device. 

1. INTRODUCTION 

In some cases, neither ray-trace analysis nor diffraction analysis can give an adequate description of an optical 
system. The designer that is faced with the problem of analyzing such a system is forced to use a ray-trace program to 
determine aberrations in the exit pupil and then introduce aberration coefficients into a diffraction model that simulate 
the propagation. We find this approach rather awkward, especially if complicated aberrations are present Our approach 
is to integrate a diffraction analysis and a ray-trace description of an optical path into one program. 

Our design is taken from a data storage application, where we must analyze the effects of objective-lens tilt The 
optical path is shown in Figure 1. A diode laser beam is collimated, and it is circularized with an anamorphic prism pair. 
The then passes through a partially polarizing beam splitter, (PPBS), which passes 70% of the irradiance polarized 
in the p direction and reflects 100% of the s polarization. The transmitted beam then reflects off of the turning minor 
into a high numerical aperture (NA) objective lens. The objective lens is usually constructed as a molded bi-aspheric 
element. The light is then focused down to the disk plane through a cover plate. At the recording layer, there is a 
complicated Interaction between the laser spot and the medium. The disk usually has a groove structure, which means 
tha t the analysis must include interaction with a periodic phase object. If the medium is a magneto-optic Sim, the 
analysis must include effects of polarization rotation upon reflection. After reflection, the beam is recollimated and 
deflected into the servo and data detection optics. 

This application requires many types of analysis. Since we ate using a laser diode, we must consider the type 
of beam tha t a laser diode generates with a physical-optics description. Interaction with the medium also requires a 
physical-optics description. Magneto-optic media require polarization effects to be simulated. A ray-trace analysis is 
necessary to quantify the aberrations introduced by the tilted objective lens. 
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Figure Z Diagram of entrance pupil to exit pupil mapping technique. 

In a previous publication Q) we have described how to analyze these effects. Our model uses a ray trace from 
the object to the image plane and a scalar diffraction propagation to simulate physical optics phenomena. There are 
commercial modeling codes, like the GLAD analysis package (2), that come close to meeting all of our requirements 
However, at the time of this writing, we are working with a code that has been developed in our laboratory. In this 
paper, we review basic principles of the model and study the effects of objective lens tilt at several planes in the optica! 
system. 


THEPRY 

Modeling of the optical system is accomplished in a two-pronged approach, as shown in Figure 2. To begin, the 
wavefront from the facet of the laser diode is propagated to the entrance pupil of the system with scalar diffraction. Next, 
we calculate a paraxial system quantifier, the ABCD matrix, which describes propagation from the entrance pupil to the 
exit pupil. Once the system ABCD matrix has been found, that information is placed into Huygens’ integral for wave 
propagation Q). To this point, we have a physical-optics description of the field in the exit pupiL 

The aberrations in the optical system must nowlie considered. We start with a point source at the laser diode 
and trace a matrix of rays to the disk plane. The optical path distance (OPD) traversed by each ray is accumulated during 
the trace. Each ray is then reverse traced to the exit pupil of the system. Ray sampling is chosen so that each ray is 
identified with a single pixel in a square matrix representing the exit pupil The square matrix of sample points is the 
same as for the diffraction propagation discussed in the previous paragraph. 

The ideal diffraction wavefront and the OPD in the exit pupil due to aberrations are combined to give the total 
exit-pupil wavefront The resulting phase is just the addition of the two component phase fronts on a pixel-by-pixel basis. 
The amplitude of the beam at the exit pupil is given by the ideal diffraction propagation of the electric field amplitude 
incident on the entrance pupil _ .... — 

The wavefront represented at the exit pupil is propagated to the disk using scalar diffraction calculations. After 
interaction with the medium, the reflected beam is propagated back through the optical system in a similar manner. 

3. MODELING RESULTS 

This model is implemented on a 486 computer that runs at 33 MHZ The modeling environment is the MATLAB 
programming language (4). MATLAB is an interactive mathematics programming environment that allows a user to 
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perform a variety of calculations, including those necessary 
for the scalar diffraction model of light propagating 
through the system. 

The setup parameters are entered through a text 
file that is edited before beginning the master subroutine. 
Input parameters for laser diode beam widths, beam 
quality, wavelength, etc, are placed in the text file. Other 
parameters, like polarization characteristics, media format, 
data patterns, etc, are also placed in the text file. The 
optical system is entered in the text file as matrix, as shown 
in Table I for the beam path from laser to disk. The input 
format for the optical system closely resembles the format 
used in other ray- trace codes. The optical system consists 
of the collimator, a beam splitter, a relay lens pair, and a 
NA=0.55 objective lens. The return path to the servo 
detectors looks very similar, except that the path is 
reversed and a lower NA lens is used to focus down onto 
the detectors. The operating wavelength for this example 
is 780 am. The stop is located at the objective lens and 
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Figure y. Sagittal (SAG) and tangential (TAN) wave 
fans for an objective lens tilted 0 JO". 
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has a diameter of 4.0 mm. 

Ideal diffraction propagation is accomplished first Next, we calculate the OPD distribution due to the ray trace. 
The field in the exit pupil is calculated from the beam emitted by the laser diode and the system ABCD matrix. A profile 
of wavefront for a lens tilted 0.5* is shown in Figure 3, which indicates a combination of coma and astigmatism in the 
sagittal (SAG) and tangential (TAN) wave fans. Next, the diffraction to the disk plane is performed. Figure 4A-C 
displays spot field distributions for 0°, 0.25°, and 0.50°, objective lens tilt, respectively. The primary aberrations are coma 
and astigmatism. The field distribution is multiplied by the appropriate disk matrices. An example of a phase matrix 
representing a groove structure on an optical disk is shown in Figure 5. Propagation back to the detector plane results 
in the irradiance distributions shown in Figures 6A-C for objective-lens tilts of 0®, 0.25°, and 0.50° , respectively. The 
small bumps located on each side of the ora&al lot>e~are due to diffraction from the groove structure. Coma, because 
it is an odd aberration, is effectively canceled when the reflected light propagates through the objective lens. Astigmatism, 
because it is an even aberration, is doubled. The irradiance distribution at the detector plane can now be summed over 
sections of the matrix to simulate effects of bi-cell, quad-cell, or other multi-element detector. 
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Figure 4. Spot irradiance distributions on the medium for objective-lens tilts of: A) 0.0°; B) 0.25°; and C) 0 . 50 °. 



4. SUMMARY 

We have shown an implementation of a combined 
ray trace and diffraction optical modelings code. Our 
example application was the determination of effects due 
to objective-lens tilt in an optical data storage device. A 
physical-optics description of the beam emitted from the 
laser diode was used, along with a physical-optics 
description of the laser interaction with the storage 
medium. Scalar diffraction theory was used to propagate 
from plane to plane. The ray-trace analysis was used to 
determine aberrations due to the tilted lens. Other effects, 
such as polarization phenomena, can also be modeled. 


Figure 5. Phase matrix representing the groove 
structure on the disk. 



Figure 6 . Irradiance distributions on the detector plane for objective lens tilts of: A) 0.0°; B) 0.25°; and C) 0 . 50 °. 
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